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Abstract 

The robustness of stellar-dynamical black hole (BH) mass measurements is illustrated using 
six galaxies that have results from independent research groups. Derived BH masses have 
remained constant to a factor of ~ 2 as spatial resolution has improved by a factor of 2 - 330, 
as velocity distributions have been measured in increasing detail, and as the analysis has 
improved from spherical, isotropic models to axisymmetric, three-integral models. This 
gives us confidence that the masses are reliable and that the galaxies do not indulge in a 
wide variety of perverse orbital structures. Another successful test is the agreement between 
a preliminary stellar-dynamical BH mass for NGC 4258 and the accurate mass provided by 
the maser disk. Constraints on BH alternatives are also improving. In M31, Hubble Space 
Telescope (HST) spectroscopy shows that the central massive dark object (MDO) is in a tiny 
cluster of blue stars embedded in the P2 nucleus of the galaxy. The MDO must have a radius 
r * 0"06. M 3 1 becomes the third galaxy in which dark clusters of brown dwarf stars or 
stellar remnants can be excluded. In our Galaxy, spectacular proper motion observations 
of almost-complete stellar orbits show that the central dark object has radius r 1 0.0006 
pc. Among BH alternatives, this excludes even neutrino balls. Therefore, measurements of 
central dark masses and the conclusion that these are BHs have both stood the test of time. 
Confidence in the BH paradigm for active galactic nuclei (AGNs) is correspondingly high. 

Compared to the radius of the BH sphere of influence, BHs are being discovered at similar 
spatial resolution with HST as in ground-based work. The reason is that HST is used to 
observe more distant galaxies. Typical BHs are detectable in the Virgo cluster, and the most 
massive ones are detectable 3-6 times farther away. Large, unbiased samples are accessible. 
As a result, HST has revolutionized the study of BH demographics. 

1.1 Introduction 

The supermassive black hole paradigm for AGNs was launched by Zel'dovich 
(1964), Salpeter (1964), and Lynden-Bell (1969, 1978), who argued that the high energy 
production efficiencies required to make quasars are provided by gravity power. Eddington- 
limited accretion suggested that BH engines have masses of 10 6 to 10 9 M Q . Confidence 
grew rapidly with the amazing progress in AGN observations and with the paradigm's suc- 
cess in weaving these results into a coherent theoretical picture. Unlike the normal course 
of scientific research, acceptance of the AGN paradigm came long before there was any 
dynamical evidence that BHs exist. 
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Table 1.1 Black Hole Mass Measurements 
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Parameters - Column 2 is the distance 



Tonry et al. 2001). Column 3 is the galaxy's velocity 



dispersion outside the sphere of influence of the BH. Column 4 is the BH mass M. with error bars 
(-^low, Afhigh) from the sources in Column 8 corrected to the adopted distance. The line with all 
columns filled in contains the adopted BH mass. Column 5 is the radius of the sphere of influence 
of the BH, r cusp = GM,ja' 2 e . Column 6 is the effective spatial resolution of the spectroscopy 
(see §1.3.1). Column 7 is the measure of spatial resolution that shows how much leverage the 
observations have on the BH detection and mass measurement. Parameters not credited are from 
Tremaine et al. (2002) or from Kormendy & Gebhardt (2001). Notes on individual objects: 

Galaxy: For Ghez (2003) and Schodel et al. (2002), a t is the pericenter orbital radius of star 
S2. Otherwise, it is the radius for the centermost radial bin of stars used in the mass analysis. 

M81 and NGC 4258: M. is adopted from Bower et al. (2000) and Herrnstein et al. (1999). 

NGC 3115: Kormendy & Richstone (1992) provide a e . The resolution a, for Kormendy et 
al. (1996a) is based on the size of the nuclear star cluster, not on the HST spectroscopy. The 
corresponding BH mass is given by the virial theorem applied to this nucleus (see their §6). 
Anders et al. (2001) modeled published data and their ground-based, integral field spectroscopy. 
Isotropic models implied M. ~ 10 9 M Q , consistent with previous results. However, they find that 
"anisotropic models reduce this to ca. 2 x 10 7 M @ ." This is inconsistent with our conclusion from 
the escape velocity argument that M. Ri 10 9 M , independent of anisotropy. Therefore, pending 
publication of the details of the the Anders et al. (2001) preliminary work, I omit this result. 

NGC 4374: I adopted M. from Bower et al. (1998), but the low-M. error bar includes the 
value suggested by Maciejewski & Binney (2001). 

For the maser galaxies, cr* is the radius of the innermost maser source used in the analysis. 
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The stellar-dynamical BH search began with two papers on M 87 by Young et al. (1978) 
and by Sargent et al. (1978). Based on the non-isothermal (cuspy) surface brightness profile 
of its core and an observed rise in velocity dispersion toward the center, they showed that 
M87 contains an M. ~ 4 x 10 9 M Q MDO if the stellar velocity distribution is isotropic. At 
about the same time, it became clear that almost no giant ellipticals like M 87 are isotropic 
(e.g., Illingworth 1977; Binney 1978) and that anisotropic models can explain the cuspy 
core and the dispersion gradient without a BH (Duncan & Wheeler 1980; Binney & Mamon 
1982; Richstone & Tremaine 1985; Dressier & Richstone 1990). Nevertheless, the Young 
and Sargent papers were seminal. They set the field in motion. 

The dynamical detection of dark objects in galaxy centers began with the discovery of an 
M. w 10 6 5 M Q mass in M32 (Tonry 1984, 1987; Dressier & Richstone 1988), a 10 7 5 M Q 
object in M 3 1 (Dressier & Richstone 1988; Kormendy 1988a), and 10 9 M objects in NGC 
4594 (Kormendy 1988b) and NGC 3115 (Kormendy & Richstone 1992). The observations 
were ground-based with resolution FWHM « 1". The BH case in our Galaxy developed 
slowly (see Genzel, Hollenbach, & Townes 1994; Kormendy & Richstone 1995 for reviews), 
for two reasons. Dust extinction made it necessary to use infrared techniques that were just 
being developed in the early 1990s. And the M, measurement in our Galaxy requires the 
study of a relatively small number of stars that are bright enough to be observed individually. 
As a result, graininess in the light and velocity distributions becomes a problem. On the other 
hand, the Galactic Center is very close, so progress in the past decade has been spectacular. 
Now the Galaxy is by far the best supermassive BH case (§ 1.3.2). 

The BH search speeded up dramatically once HST provided spatial resolution a factor of 
3 to 10 better than ground-based telescopes (see Kormendy & Gebhardt 2001 for a review). 
By now, almost all galaxies in which BHs were discovered from the ground have undergone 
several iterations of improved spatial resolution. Analysis machinery has improved just as 
dramatically. This is an opportune time to take stock of the past 15 years of progress. Are the 
detections of central dark objects reliable? Are the derived masses robust? And are the dark 
objects really BHs? The BH search is starting to look like a solved problem; assuming this, 
emphasis has shifted to demographic studies of BHs and their relation to galaxy evolution 
(see Richstone et al. 1998; Ho 1999; Kormendy & Gebhardt 2001; Richstone 2003 for 
reviews). Is this a reasonable attitude? Sanity checks are the purpose of this paper. 

1.2 The History of BH Mass Measurements 

The history of supermassive BH mass measurements is summarized in Table 1.1. 
In focusing on this history, I will be concerned with whether we achieve approximately the 
accuracies that we believe. That is, I concentrate on errors of > 0.2 dex. To what extent hard 
work can further squeeze the measurement errors is discussed by Gebhardt (2003). 

In Table 1.1, horizontal lines separate BH detections based on stellar dynamics (first 
group), ionized gas dynamics (middle), and maser dynamics (last group). All multiple 
stellar-dynamical M. estimates for the same galaxy are listed. Our Galaxy, M31, M32, 
NGC 3115, NGC 3377, and NGC 4594 have all been measured by at least two competing 
groups. M 81 has been observed independently in stars and ionized gas; both measurements 
are listed and they agree. However, consistency checks of M. values based on ionized gas 
dynamics have revealed some problems in other galaxies; these are discussed by Maciejew- 
ski & Binney (2001), Barth et al. (2001), Verdoes Kleijn et al. (2002), Barth (2003), and 
Sarzi (2003). I have not included all multiple measurements based on ionized gas dynamics. 
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1.3 How Robust Are Stellar-Dynamical BH Mass Estimates? 

1.3.1 The History of the BH Search As Seen Through Work on M32 

M32 was the first application of many improvements in spatial resolution, in 
kinematic analysis techniques, and in dynamical modeling machinery. It provides an 
excellent case study for a review of these developments. Figure 1 . 1 illustrates the remarkable 
result that BH mass estimates for M 32 have remained stable for more than 15 years while a 
variety of competing groups have improved the observations and analysis*. 

The BH in M 32 was discovered as early as possible, when the spatial resolution was 
so poor that r cusp /a* < 1. This is not surprising, given the importance of the problem. In 
astronomy as in other sciences, if you wait for a 5 a result, someone else is likely to make the 
discovery when it is still a 2 a result. The trick is to be careful enough to get the right answer 
even when the result is uncertain. Tonry (1984, 1987) got within a factor of 2.5 of the current 
best BH mass even though he made serious simplifying assumptions. His spectra did not 
resolve the intrinsic velocity dispersion gradient near the center; rotational line broadening 
accounted for the apparent dispersion gradient. Without an intrinsic dispersion gradient, his 
models were guaranteed not to be self-consistent, because there was no dynamical support 
in the axial direction. Despite this approximation, Tonry derived M. ~ (6 to 8) x 10 6 M , 
close to the modern value. Poor spatial resolution allowed considerable freedom to interpret 
dispersion gradients as unresolved rotation; since V and a contribute comparably to the 
dynamical support, trading one for the other results in no large change in M % . 

The spatial resolution of the spectroscopy improved by a factor of 30 from the discovery 
observations (Tonry 1984) to the Space Telescope Imaging Spectrograph (STIS) data from 
HST. In Column 6 of Table 1.1, the Gaussian dispersion radius of the PSF is estimated as 
follows. First, I estimate the resolution in the directions parallel and perpendicular to the 
slit as ct^ii, the sum in quadrature of the radius er* te i of the telescope PSF and of 1/2 pixel, 
and <r*±, the sum in quadrature of the radius of the telescope PSF and half of the slit width. 
The HST PSF was modeled in van der Marel, de Zeeuw, & Rix (1997b) as the sum of three 
Gaussians; for all HST observations, I use er» t ei — 0"036, the best single Gaussian dispersion 
radius that fits this sum. Finally, the effective cr» is the geometric mean of cr^y and <t*j_. I do 
not take into account slit centering errors; for some observations, these are larger than tr*. 

* The referee suggests that this result is caused by two effects that accidentally cancel because spatial resolution 
and dynamical models have improved in parallel. He suggests (1) that M. estimates increase with improving 
spatial resolution because we reach farther into the BH sphere of influence and (2) that M, estimates decrease as 
dynamical models get more sophisticated because the models have more freedom to tinker the orbital structure 
to fit the data without a BH. I disagree. (1) Reaching farther into the BH sphere of influence should not change 
M. if we model the stellar dynamics adequately well. Instead, we should get more "leverage" and smaller 
mass error bars. Of course, if we model the physics incorrectly, then more leverage may result in a systematic 
change in M. . But the change could go either way, depending on how the models err in approximating the 
true velocity anisotropy. In fact, Figure 1 .4 shows that improving the spatial resolution does not increase the 
M. values given by the Gebhardt et al. (2003) three-integral models, although it does, as expected, improve the 
error bars. For the Magorrian et al. (1998) models, improving the resolution decreases M., an effect opposite 
to that predicted by the referee. (2) Improving modeling techniques provides more degrees of freedom on the 
orbital structure, but modeling programs do not have any built-in desire to decrease the BH mass. Instead, they 
have instructions to fit the data. Again, if the real orbital structure is sufficiently well approximated by simple 
models, then making the models more complicated will not change the BH mass. And if the orbital structure is 
not well approximated by the simple models, then better models could just as easily increase M, as decrease it. 
However, the low-mass error bar on M, will decrease, for the reason the referee suggests. The high-mass error 
bar will increase. As a result, the error bars become larger and more realistic. This effect is evident in Table 1.1. 
I conclude that the consistency of M. estimates in Figures 1 . 1 and 1 .2 tells us something important, namely that 
we have been modeling the stellar dynamics of power-law galaxies well enough to derive robust BH masses. 
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Fig. 1.1. History of the stellar-dynamical BH search as seen through work on M 32: derived 
BH mass as a function of (top) publication date and (bottom) spatial resolution. Resolution 
is measured along the top axis by the Gaussian dispersion radius er* of the effective PSF 
(see text). More relevant physically (bottom axis) is the ratio of the radius of the sphere 
of influence of the BH, r cusp = GM./er 2 , to <r*. If r cusp /<7» * 1, then the measurements are 
dominated by the mass distribution of the stars rather than by the BH. If r cusp /cr* ^> 1, then 
we reach well into the part of the galaxy where velocities are dominated by the BH. Symbols 
shapes encode improvements in observations or kinematic measurements (right key) and in 
dynamical modeling techniques (left key). The data are listed in Table 1.1. 



Dressier & Richstone (1988) and Richstone, Bower, & Dressier (1990) followed with 
better observations and analysis. They fitted spherical maximum entropy models including 
velocity anisotropy. By this time, it was well known that unknown velocity dispersion 
anisotropy was the biggest uncertainty in M, measurements based on stellar dynamics. 
They were unable to explain the central kinematic gradients in M 32 without a BH. Rapid 
confirmation of Tonry's BH detection contributed to the early acceptance of this subject. 
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Since then, dynamical modeling machinery has improved remarkably. The next major 
step defined the state of the art from 1995 through 1997. This was the use of two-integral 
models that included flattening and velocity dispersion anisotropy. Essentially simultaneous 
work by van der Marel et al. (1994b), Qian et al. (1995), and Dehnen (1995) all derived 
M, = 2.1 x 10 6 M Q from van der Marel's data. Soon thereafter, Bender, Kormendy, & 
Dehnen (1996) got 3.2 x 10 6 M Q using the same machinery on CFHT data of slightly higher 
resolution. The limitation of these models, as the authors realized, was the fact that two- 
integral models are approximations. They work best for cuspy and relatively rapidly rotating 
galaxies like M 32, but they are not fully general. Still, by this time, it was routine to measure 
not just the first two moments of the line-of-sight velocity distributions (LOSVDs) — that 
is, V and a — but also the next two coefficients hj and hi, in a Gauss-Hermite expansion 
of the LOSVDs. These measure asymmetric and symmetric departures from Gaussian line 
profiles. In a transparent galaxy that rotates differentially, projection guarantees that hj, ^0. 
In general, hi is antisymmetric with V. A galaxy containing a BH is likely to have hi, > 0; 
that is, an LOSVD that is more centrally peaked than a Gaussian. The reason is that stars 
close to the BH move very rapidly and give the LOSVD broader symmetric wings than they 
would otherwise have (van der Marel 1994). Thus, as emphasized especially by van der 
Marel et al. (1994a), measuring and fitting hj, and hi, adds important new constraints both to 
the stellar distribution function and to the BH detection and mass determination. 

HST Faint Object Spectrograph (FOS) observations of M 32 were obtained by van der 
Marel et al. (1998b). These authors further "raised the bar" on BH mass measurements by 
fitting their data with three-integral dynamical models constructed using Schwarzschild's 
(1979) method. Such models now define the state of the art (see Cretton et al. 1999b; 
Gebhardt et al. 2000a, 2003; Richstone et al. 2003 for more detail). 

Finally, the most thorough data set and modeling analysis for M32 is provided by 
Verolme et al. (2002). They use the SAURON two-dimensional spectrograph to measure 
V, a, h 3 , and h 4 in the central 9" x 1 1". Also, HST STIS spectroscopy (Joseph et al. 2001) 
provides improved data near the BH. These observations fitted with three-integral models 
for the first time break the near-degeneracy between the stellar mass-to-light ratio, M/L, 
and the unknown inclination of the galaxy. Because the mass in stars is better known, the 
BH mass is more reliable. Again, the derived BH mass is similar to that given in previous 
analyses, M. =(2.9 ±0.6) xl0 6 M Q . 

So the BH mass derived for M 32 has remained almost unchanged while the observations 
and analysis have improved dramatically. It was exceedingly important to our confidence in 
the BH detection to test whether the apparent kinematic gradients near the center could be 
explained without a BH. Asked to do this, a dynamical modeling code attempts to fine-tune 
the stellar velocity dispersion anisotropy. In general, it tries to add more radial orbits near 
the center, because doing so implies less mass for the same a. Nowadays, its freedom to 
tinker is severely restricted by the need to match the full LOSVDs. However, even simple 
approximations to the dynamical structure gave essentially the correct BH mass. That is, 
M32 does not use its freedom to indulge in perverse orbit structure. The following sections 
show that this is also true in our Galaxy, M31, NGC 3117, NGC 3377, and NGC 4594. 
Dynamical mass modeling is relatively benign in such galaxies that have power-law profiles 
(for more details, see Kormendy et al. 1994; Lauer et al. 1995; Gebhardt et al. 1996; Faber 
et al. 1997; Lauer 2003). It would not be safe to assume that this result applies equally well 
to galaxies with cuspy cores. 
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1.3.2 The Best Case of a Supermassive Black Hole: Our Galaxy 

Figure 1.2 summarizes the history of BH mass measurements in galaxies with 
observations or stellar-dynamical mass analyses by different research groups. The BH case 
that has improved the most is the one in our Galaxy. Both the evidence for a central dark 
object and the arguments that this is a BH and not something less exotic like a cluster of 
dark stars are better in our Galaxy than anywhere else. 
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Fig. 1.2. Effective resolution of the best spectroscopy (top two panels) and resulting BH 
mass estimates (bottom) versus publication date. The data are listed in Table 1.1. For M 3 1 
and M32, steep rises in r cusp /<7* occur when HST was first used to observe the galaxies. 
For our Galaxy, two jumps in r cusp /<7* occur when the kinematic work switched from radial 
velocities to proper motions in the Sgr A* star cluster and when the first nearly complete 
stellar orbit in that cluster was observed. 
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A complete review of the BH search in our Galaxy is beyond the scope of this paper. Early 
work is discussed in Genzel & Townes (1987); Genzel et al. (1994); Kormendy & Richstone 
(1995), and in conference proceedings such as Backer (1987), Morris (1989), and Genzel & 
Harris (1994). Observations of our Galactic Center benefit from the fact that it is 100 times 
closer than the next nearest good BH cases, M 3 1 and M 32. For a distance of 8 kpc, the scale 
is 25"8 pc -1 . Early gas- and stellar-dynamical studies suggested the presence of a several- 
million-solar-mass dark object. In Table 1 . 1 and Figure 1 .2, 1 date the convincing case for a 
BH to Sellgren et al. (1990) and to Kent (1992). Since then, two dramatic improvements in 
spatial resolution have taken place. 

Research groups led by Reinhard Genzel and Andrea Ghez have pioneered the use of 
speckle interferometry and, more recently, adaptive optics imaging and spectroscopy to 
achieve spatial resolutions good enough to resolve a tiny cluster of stars (radius <~ 1") that 
surrounds the compact radio source Sgr A* at the Galactic Center. The Sgr A* cluster is 
so tiny that stars move fast enough to allow us to observe proper motions. This provides 
a direct measure of the velocity dispersion anisotropy. It is not large. The derived central 
mass is about 2.5 x 10 6 M©. And, even though the number density of stars is higher than we 
observe anywhere else, the volume is so small that the stellar mass is negligible. The advent 
of proper motion measurements accounts for the jump in r cusp /o-* at the start of 1997. 

A second jump in r cusp /<7* has just occurred as a result of an even more remarkable 
observational coup. As reviewed in this volume by Ghez (2003), Schodel et al. (2002), 
Ghez (2003), and Ghez et al. (2003) have independently measured several individual stellar 
orbits through pericenter passage. In the case of star S2, more than half of an orbit has 
been observed (period = 15.78 ± 0.82 years). The orbit is closed, so the controlling mass 
resides inside r peri ~ 0"0159 ~ 0.00062 pc ~ 127 AU ~ 1790 Schwarzschild radii. This 
accounts for the current jump in spatial resolution. As measurement accuracies improve, 
the observation of individual closed orbits will rapidly obsolete the complicated analysis 
of stellar distribution functions that describe ensembles of stars at larger radii. Rather, the 
analysis will acquire the much greater rigor inherent in the two-body problem. Arguably 
the orbit of S2 already contributes as much to our confidence in the BH detection as all 
stars at larger radii combined. The best-fitting BH mass, M. = (3.7 ±0.4) x 10 6 M Q , is in 
good agreement with, but slightly larger than, the value derived from the stellar-dynamical 
modeling. This leads to an important point: The above comparison in our Galaxy and 
a similar one in NGC 4258 (see the next section) are currently the only reliable external 
checks on our stellar-dynamical modeling machinery. The measurement accuracies are not 
good enough yet to show whether the models achieve the accuracies that we expect for the 
best data (±30%: Gebhardt 2003). But neither test points to modeling errors that range over 
a factor of ~ 6 as feared by Valluri, Merritt, & Emsellem (2003). 

Finally, these new observations have an implication that is actually more fundamental 
than the mass measurement. They restrict the dark mass to live inside such a small radius 
that even neutrino balls (Tsiklauri & Viollier 1998, 1999; Munyaneza, Tsiklauri, & Viollier 
1998, 1999; Munyaneza & Viollier 2002) with astrophysically allowable neutrino masses 
are excluded. The exclusion principle forces them to be too fluffy to be consistent with 
the radius constraints. Dark clusters of brown dwarf stars or stellar remnants were already 
excluded (Maoz 1995, 1998) — brown dwarfs would collide, merge, and become visible 
stars, and stellar remnants would evaporate via relaxation processes. The maximum lifetime 
of dark cluster alternatives to a BH is now a few times 10 5 yr (Schodel et al. 2002). 
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1.3.3 The Best Test of Stellar Dynamical M, Estimates: NGC 4258 

The galaxy that stands out as having the most reliable BH mass measurement 
is NGC 4258. Very Long Baseline Array measurements of its nuclear water maser disk 
reach to within 0"0047 = 0.16 pc of the BH (Miyoshi et al. 1995). The rotation curve, 
V(r) = 2180 (r/0"001)" 1/ ' 2 km s" 1 , is Keplerian to high precision. Proper motion and 
acceleration observations of the masers in front of the Seyfert nucleus are consistent with the 
radial velocity measurements along the orbital tangent points (Herrnstein et al. 1999). All 
indications are that the rotation is circular. Therefore M. = (3.9 ±0.1) x 10 7 M Q is generally 
regarded as bomb-proof. 

This provides a unique opportunity to test the three-integral dynamical modeling 
machinery used by the Nuker team (Gebhardt et al. 2000a, b; 2003; Richstone et al. 2003). 
NGC 4258 contains a normal bulge much like the one in M31 (Kormendy et al. 2003a). 
Siopis et al. (2003) have obtained HST STIS spectra and WFPC2 images of NGC 4258. 
The STIS spectroscopy has spatial resolution r cusp /(T* ~ 8.4 well within the range of the BH 
discoveries in Table 1.1. The kinematic gradients are steep, consistent with the presence of 
a BH. Three-integral models are being calculated as I write this; the preliminary result is 
that M, = (2 ± 1) x 10 7 Mq. The agreement with the maser M. is fair. The problem is the 
brightness profile, which involves more complications than in most BH galaxies. A color 
gradient near the center may be a sign of dust obscuration, and correction for the bright AGN 
(Chary et al. 2000) is nontrivial. Both problems get magnified by deprojection. 

1.3.4 A Case History of Improving Spatial Resolution: NGC 3115 

One sanity check on BH detections is that apparent kinematic gradients should get 
steeper as the spectroscopic resolution improves. We have seen this test work in M 32 and 
in our Galaxy. This section is a brief discussion of NGC 3 1 1 5 . At r cusp /er» = 59, NGC 3115 
is surpassed in spectroscopic resolution only by our Galaxy, NGC 4258, and M 3 1 . 

Exploiting the good seeing on Mauna Kea, Kormendy & Richstone (1992) found a central 
dark object of 10 9 M Q in NGC 3115 using the Canada-France-Hawaii Telescope (CFHT). 
The resolution was not marginal; r cusp /er* ~ 5.5. This is higher than the median for HST BH 
discoveries in Figure 1.3 (§ 1.3.7). Since then, there have been two iterations in improved 
spectroscopic resolution (Kormendy et al. 1996a). The apparent central velocity dispersion 
increased correspondingly: it was a = 295 ±9 km s" 1 at r cusp /cr» ~ 5.5, a = 343 ± 19 km 
s" 1 at r cusp /er» ~ 10.6 (CFHT plus Subarcsecond Imaging Spectrograph), and a = 443 ± 18 
km s -1 at r cusp /er» ~ 59 (HST FOS). These are projected velocity dispersions: they include 
the contribution of foreground and background stars that are far from the BH and so have 
relatively small velocity dispersions. However, NGC 3115 has a tiny nuclear star cluster that 
is very distinct from the rest of the bulge. It is just the sort of high-density concentration 
of stars that we always expected to find around a BH. From a practical point of view, it is 
a great convenience, because it is easy to subtract the foreground and background light as 
estimated from the spectra immediately adjacent to the nucleus. This procedure is analogous 
to sky subtraction. It provides the velocity dispersion of the nuclear cluster by itself and 
is, in effect, another way to increase the spatial resolution. The result is that the nuclear 
cluster has a velocity dispersion of a = 600 ± 37 km s" 1 . The effective spatial resolution 
of this measurement is not determined by the spectrograph but rather by the half-radius 
r/, = 0"052 ±0"010 of the nuclear cluster. This is smaller than the entrance aperture of the 
FOS. It implies that r cusp /er» ~ 59, as quoted in Table 1.1. 
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The nucleus allows us to estimate the BH mass independent of any velocity anisotropy. If 
the nucleus consisted only of old stars with the mass-to-light ratio measured for the bulge, 
then its mass would be ~ 4 x 10 7 Mq and its escape velocity would be ~ 352 km s" 1 . This 
is much smaller than the observed velocities of the stars. The nucleus would fly apart in a 
few crossing times r croS s. But r cross ~ 16,000 yr is very short. Therefore, a dark object of 
10 9 M Q must be present to confine the stars within the nucleus. 

7.3.5 A Comparison of Ground-Based and HST Studies ofNGC 3377 and NGC 4594 
Besides M 31 (§ 1.4), HST has confirmed ground-based BH detections in two more 

galaxies (Fig. 1.2). 

NGC 4594, the Sombrero galaxy, was observed with the CFHT by Kormendy (1988b), 
yielding a BH mass of M. w 10 8 7 Mq. Resolution was average for BH detections; 
7"cusp/c* = 3.8. The galaxy was reobserved with HST by Kormendy et al. (1996b) using the 
FOS atr cusp /<7* w 15.6. They confirmed the BH detection and quoted a slightly higher mass 
of 10 9 Mq. This test is weaker than those quoted above because the same research group was 
involved and because three-integral models were not constructed. However, independent 
dynamical models by Emsellem et al. (1994) agree very well with the results in Kormendy 
(1988b). 

NGC 3377 also has a CFHT BH detection; r cusp /<7* = 1.57 (Kormendy et al. 1998). The 
BH mass was M. = (2 ± 1) x 10 8 Mq. Gebhardt et al. (2003) reobserved the galaxy with 
HST at r cusp /cr* = 3.4. The improvement in resolution is smaller than normal because the 
CFHT seeing was very good and because the HST FOS aperture size was 0"2. Nevertheless, 
the improvement is substantial. Also, the analysis machinery was updated; Kormendy et 
al. (1998) fitted analytic approximations to V and a and, independently, spherical maximum 
entropy models with post-hoc flattening corrections. Gebhardt et al. (2003) fitted three- 
integral models. They obtained M, = 1-0+q 9 x 10 8 Mq, confirming the earlier result. Also, 
Cretton et al. (2003) report two-dimensional spectroscopy in the inner 6" x 3" of NGC 3377. 
Three-integral models give M. = 5.7+2 f x 10 7 Mq, corrected to our adopted distance. Again, 
the published results are consistent. 

1.3.6 Robustness of Stellar-Dynamical M. Values. I. Conclusion from §§1.3.1 - 1.3.5 

All of the ground-based, stellar-dynamical BH detections discussed in Kormendy 
& Richstone (1995) have now been confirmed at higher spatial resolution and with more 
sophisticated modeling machinery. All of the original mass estimates agree with the best 
current values to factors of 2 - 3 or better. 

Given the above tests, given the agreement between the BH parameter correlations 
implied by the dynamics of stars, of ionized gas, and of maser gas, and especially given the 
tightness of the scatter in the M, - a correlation, it seems unlikely that M. values are still 
uncertain to factors of several, as suggested by Valluri et al. (2003). Nevertheless, so much 
is at stake that we must continue to test the stellar-dynamical modeling codes. For example, 
triaxiality is not yet included. It is unlikely that triaxiality provides enough new degrees 
of freedom to greatly change the results; very triaxial configurations would have been seen 
with HST. But checking the consequences of triaxiality is under way by the SAURON team. 

All papers contain simplifying assumptions. Science is the art of getting the right answer 
using approximate analysis of imperfect data. We should not get complacent, but we appear 
to be doing reasonably well. 
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1.3.7 Application: HST BH Discoveries 

Having shown from repeat observations at better spatial resolution how well we do 
when r cusp /er* ~ 1 — 10, we now apply these results to HST BH discoveries that do not have 
repeat measurements. 
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Fig. 1.3. Spectroscopic spatial resolution for all BH discoveries in Table 1.1. Galaxies with 
repeat measurements are from Figures 1.1 and 1.2. Note that HST and ground-based BH 
discoveries have similar distributions of r cusp /cr*. However, HST has 5-10 times better 
spatial resolution in arcsec than ground-based observations (absent adaptive optics). This 
means that HST is being used to discover lower-mass BHs in more distant galaxies. 

Figure 1 .3 shows the distribution of r cusp / cr* values for all BH detections made with HST. 
It contains a number of surprises. Contrary to popular belief, HST BH discoveries are not 
being made with much better spatial resolution than ground-based BH discoveries. Only a 
few of the best HST cases have r cusp / cr* ~ 10 comparable to the ground-based BH detections 
in our Galaxy, in M 31, and in NGC 3115. On average, HST BH discoveries are being made 
at lower r cusp / cr* values than those made from the ground. Several have r cusp /cr* < 1 , similar 
to the early measurements of M 32. I am not suggesting that HST and ground-based spatial 
resolutions are similar in arcsec. HST is better by a factor of 10 (if a 0"1 slit is used) or 
at least 5 (for measurements with the 0"2 aperture or slit). What is really going on is this: 
The ground-based observations "used up" the best galaxies. For example, our Galaxy, M 3 1 , 
and M 32 are unusually close, and NGC 3115 has an unusually large BH mass fraction. So 
HST is necessarily being used on more distant galaxies or ones that have smaller BH mass 
fractions. This puts the exceedingly important contributions of HST into perspective: 

(1) HST did not find the strongest BH cases. NGC 4258 and our Galaxy were observed 
from the ground. HST observations of NGC 4258 serve to test the stellar-dynamical models. 

(2) HST has confirmed and greatly strengthened the BH cases for BH discoveries made 
from the ground. The spectroscopic resolutions r cusp /cr* for M31 and for NGC 3115 are 
now essentially as good as that for the famous maser case, NGC 4258. 
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(3) HST did not revolutionize BH detections by finding them at higher resolution. 

(4) HST has revolutionized the BH search by allowing us to find smaller BHs and ones in 
more distant galaxies. This has two important implications. 

(5a) There has always been a danger that ground-based observations would be biased in 
favor of BHs that are unusually massive. Any such bias is rapidly being diluted away. In fact, 
it was not large. Kormendy & Richstone (1995) found from ground-based observations that 
the mean ratio of BH mass to bulge mass was (M./Mb U i ge ) =0.0022+qqqq9 (they averaged 
logM./Mbuige for eight BH detections, six made with stellar dynamics and one each with 
masers and ionized gas disks). Now, the data in Table 1.1 give (M./Mb u i ge ) = 0.0013 
(Kormendy & Gebhardt 2001; Merritt & Ferrarese 2001). 

(5b) HST has made it possible to detect canonical BHs (ones within the scatter of the M. 
correlations) out to the distance of the Virgo cluster. The largest BHs can be detected several 
times farther away. This has revolutionized the subject of BH demographics. We now have 
enough detections to address the question of how BH growth is related to galaxy formation. 

(6) As Figure 1 .3 emphasizes, this subject has speeded up enormously because of HST. 

1.3.8 Caveat: Cuspy Core Galaxies 

The caveat to this rosy story is that the above tests were carried out for galaxies with 
"power-law profiles" (Lauer et al. 1995). The physical distinction between such galaxies and 
ones with cuspy cores is discussed by Kormendy et al. (1994), Lauer et al. (1995), Gebhardt 
et al. (1996), and especially Faber et al. (1997). The observations imply that cuspy core 
galaxies have more anisotropic velocity distributions than do power-law galaxies (Kormendy 
& Bender 1996). They are fundamentally more difficult for the BH search (Kormendy 1993). 
The shallower volume brightness profile p(f) gives, in projection, less luminosity weight to 
the stars in the sphere of influence of the BH. The din p/d In r term in the mass derivation 
is smaller and more easily cancelled by the effects of velocity anisotropy, which is larger 
than in power-law galaxies. Stellar dynamical BH detections in cuspy core galaxies are 
few and not well tested. Comparisons between stellar-dynamical and gas-dynamical M, 
measurements do not show universally good agreement. BH masses in cuspy core galaxies 
are more uncertain than those in power-law galaxies, and the above conclusions cannot 
confidently be applied to them. We need better tests of BH detections in core galaxies. 

1.3.9 Robustness of Stellar-Dynamical M, Values. II. What Resolution Do We Need? 
BH mass estimates made with the spatial resolution shown in Figures 1.2 - 1.3 

appear to be reliable. So how good does the spatial resolution have to be? We can now 
answer this question for two M. analysis machines, the two-integral models of Magorrian et 
al. (1998) and the three-integral models of Gebhardt et al. (2003). 

Gebhardt et al. (2003) investigate, for their objects with HST spectra and BH detections, 
how the BH mass would be affected if only the supporting ground-based observations were 
used in the modeling. The HST data are higher in resolution than the ground-based data 
by a factor of 11.2 ± 1.2. If the ground-based observations are comparable in quality to 
HST observations made with the same effective spatial resolution r cusp /<7*, then Gebhardt's 
exercise distills a clean test of the effects of spatial resolution. Modeling uncertainties are 
minimized because the same analysis machinery is used on both sets of data. Gebhardt et 
al. (2003) conclude that, when the HST data are omitted, the error bars on M, are larger but 
the systematic errors in M, are small. Here we ask how these results depend on r cusp /ov 
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Fig. 1 .4. Reliability and precision of BH masses as a function of the spatial resolution of the 
observations. The ordinate is the ratio of the BH mass as obtained from ground-based data to 
that obtained with HST kinematic data included. The error bars are from the ground-based 
data only, because I want to illustrate how estimated errors grow as resolution deteriorates. 



Figure 1 .4 shows no systematic errors in the M. values given by three-integral models, 
even at low resolution. BH masses are accurate to a factor of 1.5 or better provided that 
7"cusp/c* ~ 0.3. All BH detections in Table 1.1 satisfy this criterion. At lower resolution, M, 
can be wrong by a factor of 2 or more, but the error bars remain realistic. 

The M. measurements in Magorrian et al. (1998) have two main limitations; they are 
based on two-integral models, and they are derived from low-resolution, ground-based 
spectroscopy. They can be tested with HST spectroscopy and (mostly) three-integral models 
for 13 galaxies (open circles). When r cusp /a* * 1, the two-integral models work well; they 
underestimate the best current BH masses by a factor of 0.76 ±0.09. But when r cusp /(7» < 1, 
they overestimate the BH mass by larger factors at lower resolution. The reason for the 
systematic error is unclear. At r cusp /cr* < 0.1, M. is overestimated by a factor of ~ 5. The 
majority of the Magorrian galaxies that have not been reobserved with HST are more distant 
than the ones represented in Figure 1 .4. Therefore poor resolution plus the assumption of 
two-integral models appear to be the reasons why the ratio of BH mass to bulge mass found 
by Magorrian et al. (1998) is larger than the current value of 0.0013 by a factor of 4. 

All ground-based BH discoveries in Kormendy & Richstone (1995) had r cusp /<7* > 1 
except in the earliest papers on M32. These papers also overestimated M, (Fig. 1.1). HST 
BH discoveries made with r cusp /a* ~ 0.3 to 1 are more secure than the M 32 results derived 
at the same resolution because we now fit full LOSVDs and because three-integral models 
are more reliable than simpler models. 

Given these tests and the ones in Gebhardt (2003), it seems entirely appropriate that the 
emphasis in current work has shifted from the reliability of BH discovery to the use of BH 
demographics to study the relationship between BH growth and galaxy formation. 
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1.4 Are They Really Black Holes? 

Astrophysical arguments that the dark objects detected in galaxy centers are not 
clusters of underluminous stars (Maoz 1995, 1998) are well known. Dark clusters made of 
brown dwarf stars become luminous when the stars collide, merge, and become massive 
enough for nuclear energy generation. Clusters of stellar remnants (white dwarf stars, 
neutron stars, or stellar-mass black holes) evaporate as a result of two-body relaxation. The 
time scales for these processes are compellingly short (i.e., < 10 9 yr) only for the Milky Way 
and for NGC 4258. The next best case has been M 32 (van der Marel et al. 1998b), although 
Maoz argued that it is not conclusive. News in this subject involves our Galaxy and M 31. 

As discussed in § 1.3.2, the observation of an almost-complete, closed orbit for star S2 in 
the Sgr A* cluster restricts the central dark mass to live inside the orbit's pericenter radius, 
/"peri = 1790 Schwarzschild radii. Demise time scales for dark star clusters are now < 10 6 yr. 
Even neutrino balls are excluded (Schodel et al. 2002; Ghez 2003; Ghez et al. 2003). 

Second, M 3 1 becomes the third galaxy in which astrophysical arguments make a strong 
case against dark star clusters. Bender et al. (2003) have used the HST STIS to measure the 
velocity dispersion of the tiny cluster of blue stars (King, Stanford, & Crane 1995; Lauer et 
al. 1998; Kormendy & Bender 1999) embedded in the fainter of the two nuclei ("P2") of the 
galaxy (Lauer et al. 1993). Kormendy & Bender (1999) already suggested that the central 
dark object in M31 is embedded in this blue cluster. The STIS spectra now show that the 
velocity dispersion of the blue cluster is a = 940 ±100 km s" 1 (Fig. 1 .5). This is remarkably 
high; the red stars along the same line of sight have a velocity dispersion of only 300 to 400 
km s -1 . We can now be sure that the dark object is in the blue cluster. 

From WFPC2 photometry in Lauer et al. (1998), the half-light radius of the blue cluster 
is r/, ~ 0"06. Since all of the light of the A-type stars comes from this cluster, r;, and not the 
HST PSF or slit defines the effective spatial resolution of the spectroscopy (Table 1.1). To 
confine the stars within the blue cluster, the dark object must have a radius r, t rh- Also, M, 
is larger than we thought: the virial theorem gives M, w2x 10 8 M . This approximation is 
an overestimate if the light in the blue cluster is very centrally concentrated. However, it is 
likely that M. is at least 7 x 10 7 M . This is the value adopted in Table 1.1. 
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Fig. 1.5. Spectrum (thin line) of the central 0"2 of the blue cluster. The adjacent spectrum 
of the stars in the bulge and nucleus has been subtracted. The spectrum is continuum- 
subtracted and normalized for the Fourier correlation quotient program (Bender 1990). Flux 
is in arbitrary linear units. The blue cluster has an A-type spectrum. Heavy lines show the 
spectra of an AO V star and an AO III star broadened to the line-of-sight velocity distribution 
that best fits the cluster spectrum. This figure is from Bender et al. (2003). 
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It is important to note that M, is more uncertain in M 3 1 than in other galaxies even 
though this is the second-nearest BH case. The reason is the double nucleus. Three-integral 
models are not available and would omit important physics. Four techniques have been 
used. (1) Axisymmetric models of PI and P2 give M. ~ (4 to 8) x 10 7 M Q . (2) Models of 
the double nucleus as an eccentric disk of stars give M, ~ 7 x 10 7 M Q (Tremaine 1995). (3) 
The requirement that the center of mass of the BH and the asymmetric distribution of stars 
be at the center of the bulge gives M, ~3x 10 7 Mq (Kormendy & Bender 1999). (4) The 
virial theorem applied to the blue cluster gives M, m 2 x 10 8 Mq. These masses range over 
a factor of 7. However, all four techniques are more uncertain than three-integral models 
applied to nearly axisymmetric galaxies. An improved eccentric disk model has just become 
available; it gives M. ~ 1 x 10 8 M Q (Peiris & Tremaine 2003). The most accurate BH mass 
is likely to come from such detailed analysis of the asymmetric nucleus. Here, I adopt a BH 
mass in the middle of the above range; it should be accurate to a factor of ~ 2. 

We can now ask: Can we stuff 10 8 M Q of brown dwarfs or stellar remnants into the central 
0"06 without getting into trouble? The answer is "no" (Kormendy et al. 2003b). Following 
Maoz (1995, 1998), brown dwarfs are strongly excluded. The collision time for even the 
most massive brown dwarf (which becomes a luminous star after only one merger) is less 
than 10 9 yr. Less massive brown dwarfs collide more quickly. Dark clusters made of stellar- 
mass BHs or neutron stars evaporate in several billion years and are at least weakly excluded. 
The most viable dark cluster would be made of 0.6 M Q white dwarfs. Such a cluster would 
have an evaporation time of 10 10 yr and is not excluded by the arguments made so far. 

However, we can add a new argument. An MDO made of stellar remnants is viable only 
if its progenitor stars can safely live their lives and deliver their remnants at suitable radii. 
Progenitors get into more trouble than their remnants. They are so close together that they 
collide too quickly, as follows. The progenitor cluster must be as small as the dark cluster, 
because dynamical friction is too slow to deliver remnants from much larger radii. We get 
into less trouble with collisions if fewer progenitors are resident at one time. That is, if the 
dark cluster was made in time T, and if the progenitor star lifetime is T*, the safest strategy 
is to have T /T* successive generations, each with an equal number of progenitors. For 
T = 10 10 yr, we then calculate the time scale on which any one progenitor star collides with 
another as a function of the stellar mass and generation number. The longest time scales are 
1 8 yr for black hole and neutron star progenitors and shorter for the more troublesome white 
dwarf progenitors. Colliding stars merge and become progenitors of higher-mass remnants. 
Also, stellar mergers decrease the number of stars and increase the mass range and so shorten 
the dynamical evolution time. The result is a dark cluster with a short evaporation time. 

So astrophysically plausible alternatives to a supermassive BH are likely to fail. Our 
leverage on the M 31 BH, r cusp /a* ~ 81, is almost as good as r cusp /cr* ~ 93 for NGC 4258. 
Astrophysical arguments against BH alternatives are stronger in the latter case because we 
know in NGC 4258 but not in M 3 1 that the rotation curve is accurately Keplerian at r * <r* . 
This leads to a factor-of-ten stronger constraint on the half-mass radius of the dark object 
in NGC 4258 (Maoz 1995, 1998). Nevertheless, M31 becomes the third galaxy in which 
astrophysical arguments favor the conclusion that a dynamically detected dark object is a 
BH. This increases our confidence that all of them are BHs. 

Finally, a great variety of AGN observations, including relativistic jets and X-ray Fe Ka 
line widths as large as 1/3 of the speed of light, argue forcefully that the engines for nuclear 
activity in galaxies are BHs. 
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1.5 Conclusion 

The sanity checks that were the purpose of this paper have succeeded. Progress 
on a broad front is on the agenda for this meeting. It is gratifying to see the developing 
connection between the dynamical BH search and the AGN work that motivated it. The 
limited contact between these subjects was a complaint in Kormendy & Richstone (1995). 
Now, reverberation mapping (Blandford & McKee 1982; Netzer & Peterson 1997) has 
become a reliable tool to measure BH masses (Gebhardt et al. 2000c; Barth 2003). Ionization 
models of AGNs (Netzer 1990; Rokaki, Boisson, & Collin-Souffrin 1992) are consistent 
with other techniques (McLure & Dunlop 2001; Wandel 2002; Shields et al. 2003). The 
growing connection between BH dynamical searches, AGN physics, and the study of galaxy 
formation is a sign of the developing maturity of this subject (e.g., Kormendy 2000). The 
emphasis on BH discovery has given way to the richer field of BH astrophysics. 

Kormendy & Richstone (1995) was entitled "Inward Bound: The Search for Supermassive 
Black Holes in Galaxy Nuclei" because the BH search is an iterative process. "We make 
incremental improvements in spatial resolution, each expensive in ingenuity and money. 
[The above] paper reviews the first order of magnitude of the inward journey in radius." At 
that time, the best BH candidate, NGC 4258, was observed with a resolution cr* ~ 44,000 
Schwarzschild radii (Miyoshi et al. 1995). Now 0% ~ 23,000 Schwarzschild radii for M 31 
and for NGC 3115. The best BH case, our own Galaxy, has 0% ~ 1790 Schwarzschild radii. 
The gap between the smallest radii reached by dynamical studies and the radii studied by 
the well-developed industry on accretion disk physics is shrinking. This, too, is a sign of the 
growing maturity of the subject. 

But it is too early to declare the problem solved. Loren Eiseley (1975) wrote: 

"The universe [may be] too frighteningly queer to be understood by minds like 
ours. It's not a popular view. One is supposed to flourish Occam's razor and reduce 
hypotheses about a complex world to human proportions. Certainly I try. Mostly I 
come out feeling that whatever else the universe might be, its so-called simplicity is a 
trick. I know that we have learned a lot, but the scope is too vast for us. Every now and 
then if we look behind us, everything has changed. It isn't precisely that nature tricks 
us. We trick ourselves with our own ingenuity." 
However reassured we may be by the tests reviewed here, it is worth remembering that even 
star S2 in the Galaxy's Sgr A* cluster, which approaches to within 1790 Schwarzschild radii 
of the central engine, lives well outside the region of strong gravity. Surprises are not out of 
the question. Further tests of the BH paradigm are worthwhile to make sure that we do not 
suddenly find ourselves in an unfamiliar landscape. 
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